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Strong interaction with other particles or feedback from the medium on a Brownian particle 
entail memory effects in the effective dynamics. That motivates the extension of the fluctuation- 
dissipation theorem to nonequilibrium Langevin systems with memory. An important application 
is to the nonequilibrium modification of the Sutherland-Einstein relation between diffusion and 
mobility in the case of strong memory. Nonequilibrium corrections include the time-correlation 
between the dynamical activity and the velocity of the particle, which in turn leads to information 
about the correlations between the driving force and the particle's displacement. 



I. INTRODUCTION 

Path-integrals are robust against small perturbations 
in the dynamics and hence, make expansions easier such 
as for the derivation of response relations. That has been 
systematically applied before for the extension of the 
fluctuation-dissipation theorem to nonequilibrium sys- 
tems, at least in the Markov case [IH1]- The purpose of 
the present paper is to further extend the nonequilibrium 
linear response theory to dynamics with memory, which 
is physically often more appropriate, and which has 
been so far considered only in a "weak" non-Markovian 
case, where memory decreases exponentially fast jSJ [B] . 
Also here, the very name fluctuation- dissipation relation 
needs to be revised and perhaps altered, as the response 
obtains correlations both with the excess entropy flux 
(which is responsible for the standard relation with 
dissipation) and with the time-symmetric part of the 
action, which has been called the frenetic contribution 
as it ultimately relates to the dynamical activity in the 
process, [3]. For more practical purposes it is the latter 
frenetic contribution where the steady nonequilibrium 
forcing appears and hence, fluctuation-response relations 
can yield information about that forcing. Applications 
such as to the mobility of particles in living cells are in 
progress [7], but in the present paper we concentrate on 
the general framework, numerical exploration and some 
technical details. 

An important ingredient in the present work is the 
presence of memory in the equations of motion of a col- 
loidal particle. The origin of memory is diverse but it is 
always related to coupling with other particles and/or 
with the environment. In the case of dense colloidal 
suspensions, the reduced dynamics of a single particle 
certainly contains memory by integrating out the other 
particles. The theoretical study of the relation between 
its diffusion and its mobility is therefore advanced by 
the analysis of generalized Langevin equations (GLE). 
The latter also appear from other reduction and projec- 
tion schemes as generally treated via Zwanzig-Mori tech- 
niques [5J|S]. There, temporal scale separation or micro- 



macro transfer are the important considerations, but also 
intrinsic properties of the medium can contribute mem- 
ory effects. For the latter we have in mind visco-elastic 
media which react back on active particles from their 
previous history. These are of special interest for meso- 
scopic processes in tissues or membranes within living 
organisms which are known to respond more easily to 
external loads. 

Here we do not concentrate on the specific interactions or 
mechanism that have created the memory effects but we 
start from driven GLE for which we assume a structure 
that is relevant for a large number of cases of suspen- 
sions under the influence of external forces. See [TO] for 
a more recent microscopic-based derivation of GLE in 
the context of polymer physics. The most important el- 
ement in our modeling scheme is the principle of local 
detailed balance. It derives from the underlying microre- 
versibility which gives a strong connection between en- 
tropy flux and time-reversal breaking, [llj . As we will 
see in the next section, application of local detailed bal- 
ance leads to the so called Einstein relation, also called 
second fluctuation-dissipation relation, between friction 
and noise in the GLE |12j , even when modeling nonequi- 
librium situations. We do however not pay special at- 
tention to the choice of driving but we are interested 
here in general features and structures of the response. 
For a specific application of these general methods we 
refer to the recent work [7] on reconstructing the active 
forces from quantitative information on the violation of 
the fluctuation-dissipation theorem. We believe however 
that many more applications are waiting, in fact in all 
these cases where one can measure deviations from the 
standard Kubo-theory [T2 U3] . 

The study of response in nonequilibrium suspensions 
is of course not new, see e.g. [TIM!?] for applications 
to sheared media. In that respect the present contri- 
bution starts from GLE and investigates what are the 
general structures that determine the linear response. 
In particular the work of [TBI E2] addresses very sim- 
ilar questions and uses Martin-Siggia-Rose field theory 
to obtain a fluctuation-response theory. We concentrate 
on Sutherland-Einstein relations and we emphasize the 
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structure in terms of entropic versus frenetic contribu- 
tions, making contact with the Markov formulation in 
PQ . Close to the present work are also the results relating 
energy dissipation to the difference of the response and 
velocity correlation functions, [H [20] also for GLE. Here 
we emphasize however the modified Sutherland-Einstein 
relation connecting diffusion and mobility. 

In the next section the set-up is considered for gener- 
alized Langevin systems with Gaussian noise. They are 
driven away from equilibrium by non- conservative forces. 
We derive the linear response relations in Section |III| 
These are new results, ready to be applied in a new re- 
lation between diffusion and mobility for colloidal parti- 
cles in nonequilibrium visco-elastic media. Section [ill E| 
contains the simulation results for exploring the modified 
Sutherland-Einstein relation and adds visual information 
on the behavior of the various terms in the modified rela- 
tion. The main result of the paper is the extension of the 
work in 1] to include (even strong) memory effects and 
to be explicit also about the relevance of the correlations 
with dynamical activity and with the forcing. 



II. THEORY 
A. Set-up 

Consider the Langevin equation for the position xt and 
the velocity Vt of a (mass 1) particle in a medium at 
uniform temperature: 



dxt 
~dt 
dv t 
~dt 



(1) 



ds 7 (f - s)v s + F t (x t ) 



Vt + ht 



To lighten the notation we shall consider that, unless 
otherwise specified, integral bounds are understood to 
range from — oo to oo. We take the memory kernel 
7(f) > to be causal: 7(f) = for t < 0. The 
Markov case with friction coefficient 7 > is recovered 
whenever 7(|f|) = 27<5(f) is proportional to the Dirac 
delta function, which can be achieved for example from 
7W = 7 a cxp(— af)O(i) in the a — > 00 limit, with 
0(f) the Heaviside step function. The F t is the forcing, 
possibly time-dependent and non-conservative. It can 
include effective randomness beyond the Gaussian noise 
7]t, as e.g. in [7j. The parameter j3 is the inverse 



temperature of the environment, which we have taken 
in front of the noise rjt- The force rj t is a stationary 
Gaussian noise-process with zero mean. We wait to 
describe its time-correlations — see formula ^ below. 
The last term h t = ft 0(f) is a time-dependent (small) 
perturbation — we will linearly expand around f t = 0. 
For simplicity we use a one-dimensional notation, also 
in what follows, but the extension to other geometries 
or dimensions, sometimes essential for nonequilibrium 
effects, is straightforward. We will not use a Fokker- 
Planck description in what follows (but we use path 
integrals); actually the relation between generalized 
Langevin and generalized Fokker-Planck equations in 
the presence of position-dependent forces is not entirely 
clear and to our knowledge no results have been added 
after 1980 — see [3TJ [52] for what we do know. 

For path-space integration we need some further no- 
tation. Easiest is to take doubly-infinite paths lj — 
(x s ,v s , —00 < s < +00). The price to pay is that 
some expressions (integrals) become rather formal. We 
refer to [53] for a more detailed reference. Feasible al- 
ternatives or complements to path-integration to derive 
fluctuation-response relations for non-Markovian pro- 
cesses are known as Furutsu-Novikov theorems, see. e.g 



Because the noise rj t is a stationary Gaussian process 
the path-space measure is completely determined by the 
symmetric kernel T(t) for which 



ds r(f - s) (ri s r) r ) = 5(t ; — r) 



The weight of a path 10 is then proportional to 



(2) 



Vh{u>) oc CX P-^ / ds / dr T(r - s)rj s rj r (>!) 



rj s =v s + duj(s -u)v u - F s (x s ) - h s 



with 



Compared with the unperturbed dynamics (ht = 0) we 
have 



V h {u)=V (u)e- A ^ 



(4) 



with action 



Ah((jJ) = — — J ds J drr(r — s)h s v r — — J ds J dr J duT(r ~ s) j(r — u) v u h s 



- I ds I drT(r - s) h s F r (x r ) + 0{h 2 ) (5) 



I 

to first order in the perturbation h t - For stochastic in- tegration and path integrals in the non-Markovian case, 
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see also [33J and the more recent [57] with additional 
references. The equations define our dynamical 

ensemble for the path-space distribution with respect to 
the unperturbed dynamics ([!]). 

Finally a word about initial and boundary condi- 
tions, also important for the terminology. As should 
be clear from the start, the dynamics ([lj is not micro- 
scopic and its content depends on chosen levels of de- 
scription, including spatio-temporal scales. In the dis- 
cussion of diffusion one assumes no spatial confinement 
over the relevant time-scales. The behavior is then tran- 
sient concerning the position-degrees of freedom, but be- 
yond the inertial regime the velocities relax and become 
Maxwellian. In that regime overdamped approximations 
can be valid which can be formally obtained in what fol- 
lows by putting v — 0. For the general question of lin- 
ear response we also have in mind the case where the 
force F t = F + K t contains a time-independent conserva- 
tive force F = — V£7 from a normalizable and confining 
potential U. Under such a confinement and for time- 
independent forcing K t = K it is then assumed that 
there is a unique and smooth stationary density p(x, v) 
to which all initial data converge. We speak of an equilib- 
rium dynamics when K t — (only confining potential). 
The case of free diffusion F t — is however not strictly 
equilibrium as it needs not be stationary even in the ve- 
locity degrees of freedom. For example for free diffusion 
the Sutherland-Einstein relation will only be recovered 
when the initial data are also randomly chosen from a 
Maxwellian. In the formalism below, the important prop- 
erty of full equilibrium will be stationarity combined with 
time-reversal invariance. In fact, to the dynamics (I]) 
must still be added a relation between the noise correla- 
tions and the memory kernel, so as to ensure for example 
that for Ft = the velocities become Maxwellian. The 
next section takes this to a more general discussion by 
formulating the condition of local detailed balance. 



B. Entropy flux and local detailed balance 

The Einstein relation, also called the second 
fluctuation-dissipation relation, connects the noise cor- 
relations to the memory kernel in the friction. For an 
equilibrium dynamics, say taking F = h = in Q, 
that relation (which will follow as (10) below) is most 
simply derived from the requirement that the station- 
ary velocity distribution should be Maxwellian at inverse 
temperature /3. For our nonequilibrium dynamics (non- 
conservative forcing F) we have in general no information 
about that stationary distribution. Nevertheless, more 
fundamentally that Einstein relation arises in the weak 
coupling limit for the dynamical degrees of freedom with 
a large thermal bath as a consequence of the microscopic 
reversibility of the bath degrees of freedom. For our dy- 
namics we do keep the bath in thermal equilibrium and 
our driven particle does not react back on the bath. That 
is the basic reason why we will be able to maintain the 



standard Einstein relation. As the situation with mem- 
ory is formally more complicated and since we do not 
want to recall the details of a weak coupling analysis, 
we can rely on the so called condition of local detailed 
balance, see also [J. It amounts to assuming that, if the 
system were subjected to a confining force, it would reach 
equilibrium, and that the exchanges with the thermostat 
are the same as when the exerted forces drive it out of 
equilibrium. Formally, this means that the physical en- 
tropy flux can be recovered from the time-antisymmetric 
part of the action. That principle is indeed rooted in the 
reversibility of the underlying microscopic dynamics as 
shown in |llj . To be more specific about time-reversal 
we introduce the time-reversal operator 9 according to 
which 



0x f =x_ t , 9v t = -v_ t 
6h t = /i_ t , 0F t = F- t 



(6) 



where the last line also includes the time-reversal of the 
protocol (the time-dependence in the non-magnetic ex- 
ternal forces). When the rjt = (r)l) would be multidi- 
mensional we should also assume that the noise is time- 
reversal invariant in the sense that (^^q) = (VtVo)' -De- 
manding that the system obeys local detailed balance 
then means to require that the time-antisymmetric part 
of the action for Vq (first appearing in Q) is given by the 
entropy flux in the original (unperturbed) model; that is 
to say 

log^^(w) = -/3 Jdsv s v s +(3 JdsF s (x s )v s (7) 

The first term in the right-hand side is a temporal bound- 
ary term accounting for the kinetic energy difference be- 
tween the initial and final state of the trajectory. 

As is explicitly shown in Appendix [X] local detailed 
balance ^ is verified whenever 

(VsVt) = \ [l{t - s) + j(s -t)] = \ 7(1* - s\) (8) 

between the noise covariance and the symmetric part of 
the memory kernel. We repeat that Q is as such inde- 
pendent of F (nonequilibrium driving) or h (perturba- 
tion) as it expresses the thermal equilibrium of the bath; 
it formally appears (as Q) from requiring that the source 
of time-reversal breaking equals the (excess) entropy flux. 
The very same condition (fsb then also ensures that 



(9) 



A h (0u>) - Ah(u>) =/3 / dsv s h 



with right-hand side equal to the path-dependent ex- 
cess entropy flux towards the environment at inverse 
temperature /?, or the dissipated power by the force 
h t (setting fee = 1). The identity ([9| is also explic- 
itly discussed in Appendix [Aj we can derive ^ from 
requiring that the time-antisymmetric part of the ac- 
tion Ah is the excess entropy flux caused by the force 
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ht, thus equal to /3 J dsv s h s . We will use the notation 
S ex {uj) = A h (6w) - A h {uj) for ([9} in what follows below. 
It is worth to note that mathematically Q leads to sim- 
ply rewriting ([2]) as 



dsT(t- s)j(\s\) = 2S(t) 



(10) 



C. The time— symmetric part, or the activity 

The time-symmetric part of the action Ah is T ex (uj) 
Ah{9uj) + Ah{u) and is calculated to be 

T ex {uj) = (3 [ drH r (F r (x r ) - v r ) 



+ f J duv u J dr H r [y(u- r) - ~/(r - u)} (11) 

where we have introduced the "smeared-out" perturba- 
tion 



H r = J dsh s T(r—s) 



This time-symmetric part T ex is a function on path space 
and is also called the dynamical activity; it is related to 
the frenetic contribution in linear response playing an 
important role when away from equilibrium, see [H [3J. 
Note that the antisymmetric part of the memory kernel 
vanishes in the Markov case. In this limiting case, H s = 
h s /j and then 

TjTarkovM = ~ / ds/l * (F s (x s ) - V.) 

to linear order in ht, as before. 



III. LINEAR RESPONSE RELATIONS 

In what follows we denote by (•)/,, (•) the average over 
the paths generated by with or without the pertur- 
bation h t . The only randomness over which we average 
is the stationary noise rjt but sometimes an additional 
average over initial conditions will be mentioned. For a 
general observable local in time we write 0(x t , v t ) = O t 



A. General susceptibility 

The linear response of observable O is obtained from 
(O t )h - (O t ) = -{O t A h ) (12) 

or in terms of the generalized susceptibility \o defined 
by: 



From inserting ^ into (12) we find 



Xo(s,t) = - J dr J duT(r ~ s)-y(r - it) (O t v u ) 

+ I JdrT(r- s) ((v r O t ) - (F r {x r )O t )) (13) 

There are different ways to write that same formula. We 
can add and subtract to get the symmetric part of the 
memory kernel: 

Xo(s,t) = P(O t v a )-- I * dr [ duT(r~s)j(u~r) (O t v u ) 



J drT(r-s) {(v r O t ) - (F r (x r )O t )) (14) 



Another possibility, is to consider separately the time- 
antisymmetric and the time-symmetric part. In that case 
we follow the decomposition of the action Ah — (T ex — 
S ex )/2, and from (12 1 the susceptibility reads 



Xo(s,t) 



1 



(<TsO t 



1 



(TsO t ) 



(15) 



where 



and 



r s = _ _ 7" e 
Sh. 



Sh s 



= Pv s 



h=0 



h=0 



13 I drT{r - s)(F r - v r ) 



+ — jdu I dr v u r(r — s) ["f(u — r) — j(r — u)] (16) 



The formulation (15) separates an entropic from a fre- 



netic contribution as suggested e.g in [3j. The formulae 
(13)-(14)-(15) are the first principal resuls of the paper; 



they give a general understanding of the structure of 
nonequilibrium response also in the presence of memory 
effects. Moreover, there is the promise in each of the 
last terms in their right-hand sides to learn about the 
nonequilibrium driving exactly by the study of the 
response and especially from the frenetic contribution. 

The deviation from the Markov case is felt only in this 
excess dynamical activity T ex and not in the excess en- 
tropy flux S ex . That explains why only in nonequilibrium 
situations the fluctuation-response relations change when 
going from Markov to non-Markov; in equilibrium only 
the entropy fluxes enter in fluctuation-response relations. 
Of course, the transient diffusive case is a nonequilibrium 
situation, and we should be careful when possibly iden- 
tifying F r — with the equilibrium case. 



B. Consequence of causality 



h=0 



Causality requires that observations before a certain 
time are not influenced by perturbations after that time. 
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As a consequence, from (15), 



(a u O r ) = (r u O r 



(17) 



C. Equilibrium dynamics 

1. Confined case 



for a time-ordering u > r. But suppose now that the 
averages satisfy time-reversal invariance so that 



The equilibrium limiting case can be achieved when- 
ever the force field F derives from a potential function. In 
that case time-reversal invariance applies, and one should 



{r s O t ) = sgn O (r_ s O_ t ) , sgn O (<7_ s 0_ t ) = - (a s O t ) recover from (fl2j) the standard fluctuation-dissipation 



where sgn O is the parity of observable O under time- 
reversal. Then, under that time-reversibility and as a 



result of the causality relation (17) for u = — s > r = —t, 



(a s O t ) - (r s O t ) = (a.O t ) - sgn O (r_ s O_ t ) 

= (a s O t ) - sgn O (<7_.0_ t > (18) 
= (a s O t ) + (a s O t ) = 2/3 (v s O t ) 



which, upon inserting in (15), yields the standard 
fluctuation-dissipation relation 



XO (s,t)=P(O t v s )@(t-s) 



(19) 



The next section comes back to this with yet another 
derivation. 

Another consequence of causality is that the action |5| 
verifies 



(O t A h 







theorem (19). To check this more explicitly, we first con- 



sider the response of an observable O t = 0(x t ,v t ) which 
is even under time-reversal (i.e. OOt = O-t). In this 
case 



(O t v^ u ) = - (0- t v u ) 

(F r {x r )0- t ) = {F_ r {X- r ) O t ) 
(v r O-t) = (i-r O t ) 



(22) 



In Appendix [B] we show that the time anti-symmetric 
part of the response \ yields 



Xo(s, t) - Xo(t, s) = & {O t v s ) 



(23) 



By causality, this automatically implies the fluctuation- 
dissipation relation (19). 



Likewise, when the observable Ot is odd under time 
reversal (i.e. OOt = —O-t), the right-hand sides in 
Eqs.(22 ) have to be multiplied by —1, and thence one has 



to consider the time-symmetric part Xo{s, t) + Xo(t,s) 
in order to indeed recover (19). 



when h s — for s < t. That immediately implies that 
for all s > t, 



2. Free diffusion 



dr J duT(r — s) j(r — u) (OtV u ) 

= [ drT{r- s) {(F r (x r )O t ) - (v r O t )) (20) 



There is a simpler identity that applies when the original 
dynamics is time-homogeneous, i.e., when F t = F does 
not explicitly depend on time t. Then we can think of 
the unperturbed averages (•) as a steady regime. In that 
case we take time t very negative in ( 20 ) , multiply both 



sides with (rj s r] w ) for arbitrary w > t, integrate over all 
s and use the identity ([2| to obtain 



+oo 



du^(w - u) {O t v u ) = (F(x w )O t ) - {v w Ot) 



In the Markov case this identity is, for w > t: 



l^(O t x w ) = (F(x w )O t )~(v w Ot) (21) 

which is readily recognized as (O t T) w ) = for white noise 
r] w ,w > t. 



We next consider the case where F r = and there is no 
confinement on the relevant time-scales. In this case the 
velocity of the particle relaxes alright to a Maxwellian 
steady-state in the long time limit, but its position dif- 
fuses. (It could be anomalous diffusion for slowly decay- 
ing kernels.) Therefore, starting from a fixed position and 
Maxwellian velocity, the velocity response Xv satisfies 



Xv(s,t) =Xv(0,t~s) = /3(v s v t ) 



(24) 



for s < t, but the position dynamics remains in the tran- 
sient regime. It is however possible to recover a formula 
similar to ( |19[ ) which involves the mean square displace- 
ment Ax 2 (t) — ((x s +t ~ x s) 2 ) as we now explain. 
We consider the case where the observable O is the po- 
sition x. The position response 



Xx(t) = ^{ x t)h 



h=0 Jo 



h=0 



satisfies 

d_ 
d7 



t p-t 

Xx (t)= I ds X v(s,t)=f3 / ds (v s v t ) (25) 
o Jo 
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On the other hand, 



From the general linear response (15 1 we know that 



Ax 2 (t) 



s+t ps+t 

du / dr (v u v r ) 



d f s+t 



As a result, we get the equilibrium-like result 



(26) 



which, for the Markov case, is an identity attributed to 
Virasoro in [55]. The fact that this relation remains sat- 
isfied also in the presence of strong memory is directly 
caused by F = and that only the entropic contribution 
matters. 



D. Modified Sutherland-Einstein relation 

We now show how to connect the mobility of the parti- 
cle (which is related to the velocity response to a constant 
force), and the diffusion properties, related to the time 
behaviour of the mean squared displacement. Remember 
that the perturbing field is a step function h t = h®(t), 
with constant field h. The time-dependent mobility is 
then defined as: 



M(t) 



1 d_ 
t dh 



((x t -x )) h 



(27) 



h=0 



M(t) 



JO 



ds / drxv{s,r) 



(28) 



The time-dependent diffusion coefficient is defined as 



D(t) = ^Ax 2 (t) = ^ I dr I d.s /,.., , v > (2<>) 



(v s v t ) 



2 , s 1 , 

-gXv{S,t) + -g {T S V t ) 



(30) 



which can be replaced in ( 28 )-( 29 ) to give a relation be- 
tween M, D and r, 



M(t) = /3D(t) - 



t r-t 
dr j ds (t s iv) 



We see how the violation of the Sutherland-Einstein rela- 
tion M = (3D is related to the time-averaged correlation 
between displacement and dynamical activity, 



M(t) = (3D{t) 



1 

2f 
1 



ds (r s (x t - x )} 



(31) 



0D(t) -^{{xt- x ) 



This formula has a general validity, independent of ini- 
tial conditions. The time-averaged dynamical activity 
due to the perturbation has become here an important 
observable for describing the deviation from the stan- 
dard Sutherland-Einstein relation. Since t s is time- 
symmetric, the last correction term will vanish when 
time-reversal symmetry gets established; that happens 
in the case of 

free diffusion upon averaging over the initial equilib- 
rium distribution for velocities (as in the previous sec- 
tion). In general however we can further detail the ex- 
pression by using the explicit form of t s from ( 16 1, yield- 
ing three terms in the correction: 



M{t) = f3D(t) + d(t) + C 2 (t) + C 3 (t) 



(32) 



As we are however most interested in the nonequilibrium 
situation with possible drift, it is useful to discard this 
effect by considering a slighly different definition of the 
diffusion coefficient, relabeled here as D: 



x ); {x t - x )). 



(33) 



The notation (A; B) refers to the connected (also called 
truncated) correlation function (AB) — (A) (B) . The ana- 



logue of the modified Sutherland-Einstein relation (32) 
is explicitly given by: 



M(t) = PD(t) + g J ds j drT(r-s)((x t -x ); v r ) 



Ci(t) 



4f 



(3 r 



+ 77 / ds U dr / dwr(r - s){7(r - u) - >y(u - r)} ((x t - x );v u )] -— / / drT(r - s) (F r (x r ); (x t - x )} ds 



2* Jo 



c 2 (t) 



c 3 (t) 



(34) 
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This is the third key relation of the present paper (after 
the general ( 13 1— ( 14 )— ( 15 1 and the equilibrium-like (p6|). 



We have postponed a detailed derivation of this result 
in Appendix O Notice that the functions Ci of Eq.(32) 
are simply given by the functions Ci with a standard 
correlation function instead of a connected one. More 
concretely, if d = (A; B) then d = (AB). We shall 
further explore this relation numerically in the next sec- 
tion. The terms C\ and C3 are also present (although 
without convolution with T) in the Markov case, and ac- 
counts for nonequilibrium effects due respectively to the 
inertia and to the nonequilibrium forcing. The term C\ 
vanishes outside the inertial regime. The term C3 is most 
important for the inverse problem of reconstructing the 
nonequilibrium driving from violation of the Sutherland- 
Einstein relation. The term Ci has no analogue in the 
Markovian case, and is hence only due to memory effects 
under nonequilibrium dynamics. 



force F. Again we checked in all cases that for exponen- 
tially decaying memory kernels "/(t) the results of [T] are 
reproduced. We concentrate on power law decay and we 
consider finite times. In Fig|3]we consider a vector force 
on the plane to induce vortices, similar to the rotational 
force in [I]. We take F = Ag, with A the amplitude, 
where 



g x (x,y) = a{r - y/2)(y - i) 
g y (x,y) =a{r-V2)()- - x) 



(35) 



a — (1 — 2^2,a:mod3)(l — 2(5i y moc j 2) 



E. Examples 



We present here simulation results for the dynamics 
for some three choices of the driving F. The method 
to generate colored noise is outlined in Appendix [D] For 
the evaluation of the terms in ( [34] ) we note that T ap- 
pears as a convolution and the Fourier transform of T is 
computed from its definition ([2]). 

The purpose is to visualize the various terms in the mod- 



ified Sutherland-Einstein relation ( 34 1 , similarly to the 



work in [T] but with the extra ingredient of memory. Of 
course there is no strong need to verify or to confirm the 
mathematical formula ( 34 ) as the linear expansion is suf- 



ficiently controlled, but it is somewhat interesting to see 
the contributions of the various terms in ( 34 ) and to see 
what influences them. 

We start with the free diffusion (F — 0). The sys- 
tem is diffusive if, for large times, D(t) reaches a limit 
(the diffusion constant) D. On the other hand, when 
the memory gets long time tails, with j(t) decaying al- 
gebraically, the diffusion can become anomalous, [2"5] . 
We take two different examples, j(t) = 1/(1 + t) and 
7(i) = 1/vl + 1. As results show, both the diffusion 
and the mobility follow a corresponding temporal behav- 
ior, with for the first example D(t) ~ 1/logt and for 
the second example D(t) ~ (subdiffusive motion). 

The plots in Fig[l] and Fig(2] show the subdiffusive be- 
haviour due to the memory effect. More generally for 
free diffusion D cx (L j(s) ds) _1 for large t. Note that 
we treated free diffusion with some fixed initial condition 
(without initial velocity-averaging) so that the standard 
Sutherland-Einstein relation gets established only in the 
long time limit. Note that the inertial regime is rather 
short-lived, C\ rapidly being very small, but the memory 
effect as present in Ci postpones the equality between 
f3D(t) and M(t) to longer times. The sign of Ci carries 
no special information. 

For the other examples we switched on some rotational 



for distance r = J~{x — ^) 2 + (y — \) 2 . The somewhat 

involved definitions assure that the particle does not 
undergo a net drift; the forces are purely rotational 
and not translational now. To make sure, mobility and 
diffusion are now matrices but the off-diagonal elements 
are approximately zero. We also find that the diffusion 
in the x-direction is bigger than in the y-direction. 
Moreover, for bigger A, the diffusion increases, while the 
mobility remains almost constant (and even somewhat 
decreases). 

Finally, we consider the result of driving in 1 dimen- 
sion. The nonequilibrium force is obtained from a 
periodic potential, which is like confining the particle 
to a toroidal trap, and adding a constant field. In 
formulae, the nonequilibrium force is F(x) = A + sinx, 
where A is a constant. Fig |4] and Fig (5] show the result 
again for the two long-memory kernels. As we have 
found in the previous section the relation between 
diffusion and mobility gets modified. Our simulations 
confirm in all cases that the diffusion depends on the 
external forcing more strongly than does the mobility. 
In the power law decaying memory for A = 1.5 the 
increase of the diffusion is seen. C\ is still almost 
zero, Ci vanishes in the long time and the diffusion 
and mobility are not proportional any more; there 
is now also the essential term C3, which is negative 
because of the positive correlation between force and 
displacement. However, for stronger memory and with 
A = 1.5 there is little difference with the case / = 0. 
To make the contribution of C3 more prominent we 
have taken a larger driving force, like F(x) = 8 + 8 sin a; 
as in Fig[6j we see that C\ is still zero, Ci is getting 
smaller faster than before and C3 is more important now. 




Time Time 



FIG. 1. The time-dependence of the mob ility M, the diffusion 
D and the corrections Ci and C2 of ( 34 1 after multiplying all 
with logt. Here we have long memory 7(2) = 
diffusion F = at = 1. 



Y^-j: and free 




1500 2000 
Time 



FIG. 2. Same case of free diffusion as in Fig. [T] but with a 
time-rescaling of ^ft for memory kernel j(t) — ^f+i - Still 
F = 0,/3 = l. 



FIG. 3. Diffusion and mobility in the case of a two- 
dimensional rotational forcing (351 for memory j(t) = 
and amplitudes A = 10, A — 20. 




500 1000 1500 2000 2500 
Time 



FIG. 4. Same set-up as in Fig. [T] but for a forcing F(x) = 
1.5 + sin a; (external field over periodic potential) at inverse 
temperature j3 = 1. The rescaling of the mobility, the diffu- 
sion and the corrections in ( 34 1 is by multiplying all with log t 
for memory kernel j(t) = 



T 

1+t 



IV. RELEVANCE AND CONCLUSION 

The fluctuation-dissipation theorem has an extension 
to nonequilibrium generalized Langevin systems which 
preserves the splitting of the response in an entropic 
and a frenetic contribution. That general statement 
is of course relevant in today's search for useful linear 
response formula away from equilibrium. The presence 
of memory is especially relevant for dense colloidal 
suspensions. We have here considered a driving which 
can be time-inhomogeneous or even random but we 
have also assumed the presence of Gaussian (correlated) 



noise to start an expansion from path-integrals for an 
underdamped dynamics. The Gaussian correlations are 
in fact connected with the memory kernel in the friction 
via the condition of local detailed balance. An important 
application is to the extension of the Sutherland-Einstein 
relation between diffusion and mobility. Because of the 
nonequilibrium condition, the diffusion constant is no 
longer alone in determining the transport properties 
of colloidal particles. That was in particular seen in 
detailed simulations for various nonequilibrium diffu- 
sions, possibly anomalous, in particular for exploring 
the role of the nonequilibrium forcing and the influence 
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of memory. An interesting conclusion is that the 
nonequilibrium corrections to the Sutherland-Einstein 
relation are related to the time-correlations between 
the so called dynamical activity and the velocity of the 
particle, which in turn leads to information about the 
correlations between the driving force and the particle's 
displacement, [7j- Such an analysis provides a more 
general framework for discussions on the violation of 
the fluctuation-dissipation theorem and is an alternative 
for the use of effective temperatures. We believe that 
both the direct question (predicting the response) 
as the inverse question (deriving information on the 
nonequilibirum forcing) can be attacked within the given 
formalism. 



500 1000 1500 2000 2500 3000 
Time 



FIG. 5. Same set-up as in Fig. [2] but for a forcing F(x) = 
1.5 + sin a; (external field over periodic potential) at inverse 
temperature j3 = 1. The rescaling of the mobility, the diffu- 
sion and the corrections in ( 34 1 is by multiplying all with vt 
for memory kernel j(t) — ^7j=p- 
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FIG. 6. The rescaling of the mobility, diffusion and the 
corrections by multiplying with \[t for "/(t) = ^y^+t 1 an< ^ 
F(x) = 8 + 8sinx,/3 = 1. 



We start with deriving ([7j), From writing out the action 
log gp^(^) in the same way as for ^ we find that 



dVov 



ds J drT{r — s)v r J du py(r — u) + 7(1* — r)] v u 
dsjdrTtr-s) jdu l[s -u) Vu jd Wl{ r- W ) Vw 
ds J drT(r — s) J du~/(u — s)v u J dwf(w — r)v w 
ds / dr r(r — s) F r (x r ) / du [y(r — u) + j(u — r)] v u 



(Al) 
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which can be rewritten as 



ds J drT(r—s)v r J du[j(r — u) + "f(u — r)]v u 
ds J drT{r ~ s) J du^/(s — u) v u J dw~/(r — w)v u 
ds J drY(r — s) J dwy(u — s) v u J dw~f(r — w)v u 
ds J drT(r—s) J dwf(u — s) v u J dwj(r — w)v n 
ds / drT(r— s) / du^(u — s)v u I dwj(w — r)v ll 



(A2) 



2 

from which we arrive at 



+ — I ds I drT(r — s) F r (x r ) I du [y(r — u) + j(u — r)] v u 



ds / drT(r — s) du dw [j(s — u) + j(u — s)]-f(r — w) v w v u 



log Jj^M = ~% I I iIr T{r - - ' ''<■ ! <>" :. { r - 11 1 + ~ 1 " -- r] \ <•« 

_ P 

4 

+ — j <ls j dr F( /• — ••< ) / (III j ihr [~ ( /• — ir ) -I- " ; ( (/• - /')]-.(// ■ - ! r, r c„ 

I 
2 



+ — I ds I drT{r — s) F r (x r ) I du[j(r — u) + j(u — r)]v u (A3) 



Since (r] s r]t) = \ [j(t — s) + j(s — t)], as implied by local finally recovers Q. 

detailed balance ([8]), and by using the definition of the Likewise, when we add a time-dependent perturbation ht 

symmetric kernel T(t), the second term and the third we have 
term of the right-hand side cancel each other and one 



ds / drT(r — s) I du I dw [7(5 — u)j(r — w) — j(u — s)j(w — r)]v u 



log jjj^M = — g / '/>• / '//• Hr M r, / iln [Mr - n) + ~Mi - r) r„ 

_ P 
4 

+ — I lis I 1I1 Tir ■ s) F, ( ,r, ) I <lii [~[r •</)••• Mm- r) r„ 

P 
2 

(A4) 

The excess entropy flux is then 

S°*(u) = A h (6u>) A h (u) = log - log 



+ — I ds j drT(r — s)h s j dw["f(w — r) + j(r — w)] v v 



ds / dr r(r — s) h s / dw[7(w — r) + j(r — w)] v v 



which indeed ensures ^ upon using (10). Appendix B: Recovery of the equilibrium 

fluctuation— dissipation theorem 



We show how to get from ( 12 ) to ( 23 ) when ( 22 1 holds. 
We start by writing 
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r and u' — u one gets 



Xo(s,t)-xo{t,s) = - / dr / d«r(r-s)7(r-w)(O t v„) 
+ drT(r - s) {(v r O t ) - (F r O t )) 



drT(r-t)((v r O s )-(F r O s }) (Bl) 



Xo(s,t) - xo(t,s) 



dr J duT(r — s)"f(r — u) (OtV u ) 
dr / duT(r - s)j(u ~ r) (O s v t+s - u ) 



J drT(r - s) [(v r O t ) - (F r O t ) 

- (Vt + s-r0 s ) + (F t+s _ r O s )} (B2) 



Finally, by using the time-translation invariance of cor- 
relation functions plus the time-reversal conditions ( 22 ) 
one will observe that the square brackets term in the last 
integral vanish and the first two integrals simplify thanks 
to ([8]) and yields to Eq.Q. 



In the last two integrals we can perform the change of 
variable r' = t + s — r and v! = t + s — u. By relabeling 



Appendix C: The modified Sutherland-Einstein 
relation in connected form 



By inserting the dynamical activity in (31 1 we obtain 



*7, 

f) " 



ds / dr r(r — s) ((xt — Xq) v r ) 



+ 4t J ds I dr I dur ( r - s ){"/( r -u)- 7(u - r)} ((x t - x ) v u ) 



2t 



drT(r — s) (F r (x r ) (x t — xq)} ds 



We now rewrite this relation as follows 
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M(t) = ^.{(xt - x ) 2 ) - ~{x t - x ) (x t ~ x ) + ^,( x t- x Q ) (x t - Xq) 



+ 



ds / drT(r — s) ((xt — Xq) v r ) 



ds / drT(r — s) (xt — Xo) (v r ) 



I 
2* Jo 



ds / dr r(r — s) (xt — x§) (v r ) 



s){7(r — u) 




r)}((xt- 


- x ) v u ) 


s){j(r — u) 




r)}(x t - 


xq) (v u ) 


s ){l{ r ~ u ) 


-l(u- 


r)}(x t - 


x ) (v u ) 



drT(r — s) (F r (x r ) (x t — Xq)) ds + — f I drT(r — s) (F r (x r )) (x t — Xq) ds 



2t Jo 



drT(r — s) (F r (x r )} (xt — xq) ds 



(CI) 



Following the definition of the connected correlation 
function we arrive at 



M(t) = (3D(t) + C\(t) + C 2 (t) + C 3 (t) 

+ •K,( x t ~ x ) (x t - Xq) + — / ds I drT(r - s) (x t - x ) (v r ) 



2t It J 



+ if ' ^ ' ^ ' ~ S )i^( r ~ U ) ~ "f( u - r )} (Xt - Xq) (v u ) 



I 
K Jo 



drT(r — s) (F r (x r )) (x t — xq) ds 

(C2) 



r 



From the Langevin equation we have 

du-f(r - u) (v u ) + (F r (x r )) 



cellation for the forcing term is clear; and the other terms 
follow as: 



After substituting this in (C2), all the terms cancel out 
each other and only the first ine will remain. The can- 
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— J ds J dr J duT(r — s) j(r — u) (v u ) = — — J ds j dr j duT(r ~ s) j(r — u) (v u ) 



t 
4 


4 

t 
4 Jo 



ds dr duT(r — s) j(r — u) {v u ) 



At 



ds J dr J duT(r — s) j(u — r) (v u ) 
ds J dr J duT(r — s) j(u — r) (v u ) 

ds dr dvT(r — s){"f(r — u) + 7(1* — r)} (v u ) 



4* 7o 



ds dr duT(r — s){7(r — u) — 7(1* — r)} (v u ) 



-(x t - Xq) 



I 
At 



ds j dr j duT(r — s){7(r — u) — 7(1* — r)} (i> u ) 

(C3) 



Appendix D: Simulation of colored Gaussian noise 



where A is the sampling interval of time. Finally, is a 
Gaussian complex random number in Fourier space with 
zero mean and correlation, 133 , 



We sketch the algorithm to numerically generate a 
stationary Gaussian coloured noise £t for a given time- 
correlation function 7. The strategy is to transform to 

Fourier space as e.g in [30JI3I], then to simulate it there Tq generate a Gaussian com plex random number with 



S. 



fj,,N-vi a fi — a N-/i 



(D4) 



as e.g. in [32], and then to transform the solution back to 
real space. We next explain that scheme in more detail 
Suppose the noise satisfies 



<*(*)> =0 

mm)=~f(\t-t'\) 



(Dl) 



In discrete Fourier space the noise can be constructed as 

i(u M ) = y / iV 7 K) a M , M = 0,...,iV-l (D2) 

with N even, and £(w M ) and 7(w^) being the Fourier 
transforms of and 7(f), respectively; the w M is defined 
as 



lu u = 2ir— — — 
A NA 



UJ N . 



correlation given in D4 1 , we write a M = a„ + ib^ in terms 
of its real and imaginary parts: = + ib^ if /1 > N/2 
and otherwise = aN-^ — ibN-n- Here, a and b are two 
Gaussian real random numbers which are uncorrelated 
and have zero mean and covariance 



(buK) = - (S^u - S^N-v) 

(D5) 

It is then straightforward to do the inverse Fourier trans- 
form 



N-l 



^=0 



(D3) and to see that the correlations reproduce (Dl I. 
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